
NPS ARCHIVE 
1958 

ANDERSON, L. 



ATTENUATION OF REPEATED 
SHOCK WAVES IN TUBES 



LAWRENCE R. ANDERSON 
AND 

DONALD J. MERMENS 






LIWARV 

US. NAVAU POSTGRA'.i'J.M* 
mont:«ey, calif ' 



C T 



* 






* 



I 



t 

ATTENUATION 

of 

REPEATED SHOCK HAVES 

in 

TUBES 



* 



by 



Lawrence R. Anderson 
Lieutenant Colonel, United States Army 



and 



Donald J. Me hr tens 
Major, United States Army 



Submitted in partial fulfillment of 
the requirements for the degree of 

MASTER OP SCIENCE 
IN 

PHYSICS 



United States Naval Postgraduate School 
Monterey, California 
1958 



iMrs /- \c w vc 
YVb<J> 

I t- ; y > Y c\: ■ ■ 



I3KARY 

• A1 °C iCRADUA ,-jql 
•-ON' 1 ^M'fOP .A 



ATTENUATION 

of 

REPEATED SHOCK WAVES 
in 
TUBES 

by 

Lawrence R. Anderson 1 
and 

Donald J. Mehrtens 
This work is accepted as fulfilling 
the thesis requirements for the degree of 
MASTER OF SCIENCE 
IN 

FflYSICS 
from the 



United States Naval Postgraduate School 



ABSTRACT 



The field of attenuation of repeated shock waves bears 
much investigation. In this field there are theories to 
explain the attenuation of infinite plane waves and infini- 
tesimal waves but no suitable theoretical explanation of the 
attenuation of finite repeated waves, Professor I. Rudnick 
has proposed an eqtiation to explain these phenomena but has 
failed to obtain experimental proof of its validity. Some 
work has been done in this field by other investigators. 

This investigation was under taken to determine the validity 
of the theory over a series of tubes of various diameters 
and at several frequencies. 

This investigation has failed to substantiate the 
validity of the equation but has instead produced results 
which support the previous work done in the field which 
tend to negate the theory. 

The investigators wish to express their than’rs to 
Professors II. I'edwin and C. B. Wilson, Jr. of the U. S. 

Haval Postgraduate School for their encouragement and help- 
ful suggesticns during the conduct of this investigation. 
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1. Introduction 



Shock waves in gases are characterized by a pressure 
discontinuity at the shock front. In the field of gas 
dynamics, three distinct types of shock waves are observed: 

(1) Standing shock waves such as are generated from the 
leading edge of an aircraft flying at supersonic speeds, 

(2) single traveling shock waves which are generated by an 
almost instantaneous change in pressure differences such as 
occur in an atomic explosion, and (3) repeated traveling 
shock waves which are generated by intense sound waves. 

Much investigation has been undertaken for standing shock 
waves and the theory is rather complete. The Rankine- 
Hugoniot relations sufficiently describe the single traveling 
shock wave at the present time. Theory exists for the in- 
finite plane repeated shock waves but an adequate theory 
does not exist for finite plane repeated shock waves. This 
investigation is confined to the attenuation of finite plane 
repeated shock waves in tubes. 

Shock wave research is important to ascertain basic 
data and methods available for dissipating exceptionally 
high noise levels without incurring large back pressure to 
engine exhaust, etc. This is of particular importance in 
this jet age. 
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2. Characteristics of repeated shock waves. 

Both single shock waves and repeated shock waves can 
be generated under controlled laboratory conditions. The 
usual method for generating a single shock wave is by the 
rupture of a diaphram separating gases at two different 
pressures. This shock wave propogates with a velocity 
greater than sound and is the familiar type generated by 
explosions, atomic blasts, etc. Figure 1 depicts this type 



of shock wave. 





Fig. #1 

Repeated shock waves travel with the normal speed of 
sound and have a saw-tooth wave form which is stable ex- 
cept for a gradual amplitude attenuation due to absorption 
processes. The shock front is formed by the distortion of 

i 

a high amplitude progressive sound wave (whose initial 
pressure variation is in the form of a sine function) due 
to the peak overtaking the trough. As the crest of the sound 
wave travels at a velocity equal to the sum of the particle 
velocity and the velocity of sound and the trough travels at 

^G. IVerth, "Attenuation of Repeated Shock ;/aves in Tubes", 
UCLA 19f>3 
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Pressure 



a velocity less than the velocity of sound, the gradient in- 
creases on the leading front of each wave and decreases on 
the trailing front and the shock is formed. See Pigs. 2a, 
b 2c c below. 




Pig 2a High Amplitude Sound Wave 





Pig 2b Shock Wave Variation Pig 2c Shock Wave Variation 
of Pressure with of Pressure with Time 

Distance 
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3. Attenuation of reperted shock waves theory. 

Professor I. Rudnick of UCLA, in his Technical Report 
to the U. S. Navy [lJ has developed the best formula to date 
explaining the attenuation of repeated shock waves. By 
using the Rankine-Hugoniot relations and e ousting the in- 
crease in^ entropy across a shock front with the rate of 
energy decrease, he derived the following formula: 

^ ~i •' (*-*•) 



vjhere 



6 



= P 2 ~ Pi 
Pi 



A = wave length 

X - distance from 
sound source 



r= C p /C v (ratio of Specific Heats) 

The above formula assumes a saw tooth wave form and is 
one dimensional implying the infinite plane wave. The 
formula also ignores the effect of tube boundaries. No 
adequate theoretical treatment for finite amplitude or re- 
peated shock wave s in tubes exists. The be3t available for 
the former is Kirchoff’s infinitesimal solution,^ where 

i 4,(s)' ('W'Pfa ' r -r (&*]*(?)* 

and /< - coefficient of viscosity 
f = frequency 



A = thermal conductivity 
P = radius of tube 



^J. Acoustical Soc of America, Vol 26 , No 1, Jan 1951(., 
p. 61 
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V 




Wilson and Bies in a similar report to the Navy [2j 
used Rudnick's basic equation which combined the tube pnd 
the infinite plane solutions and accounted for the looses 

in the formula : 

= a-j_x * 



in the tube walls^ 

I 



* - 14 

where l/r 

' 0 A 



v 's a 2 x 

ave. 



a l 

a 2 

c 

X 

r 



id/r - i/a ) 

1 + r (attenuation at shock front) 



Zt 

/ 7 0 Sfv ‘ 
d c 



(additional losses at the 
tube walls ) 



= velocity of sound 

« distance from sound source 

= Cp/C v (ratio of specific herts 
which for air is 1 . 14 . ) 



A = ’Wave length 
d = diameter of pipe 
f = frequency 
P2 = peak pressure 
Pi » trough pressure 

The quantity v' is a function of the kinematic coefficient of 
viscosity and the thermometric conductivity coefficient for 
air.^ Wilson and Bies [ 2j found that the term involving a.£ 
is at most only 5$ of the term involving aq for their experi- 
ment. 



-^Wilson and Bies Technical Notes Report 79» Soundrive 
Engine Company, Los Angeles, California, 1953* Pg 

^■Rayleigh, ’'Theory of Sound”, Dover Publications, 19l|-5> , 
Vol II, pp 3?3-325 
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4. Previous Investigations 

Previous investigations of the attenuation of shock 
waves in tubes have been carried out in three reported 
cases. The first such investigation was by Rudnick and 
Leonard under a contract by Soundrive Engine Company with 
the Office of Naval Research and reported upon in 1953 



The investigation involved the attenuation of extremely 
high amplitude sound pressures generated by a siren with the 
air supply furnished by utilizing two Packard- Kerlin aircraft 
engine superchargers. Rudnick and Leonard used a 60 foot 
long pipe with a ten inch inside diameter for the investi- 
gation and studied frequencies of 20-200 cycles per second. 
Using the formula 



a graph was plotted of 1/^- vs x. Theoretically such a plot 
should yield a straight line of slope_£lL* However, the ob- 



served data when plotted in this way had slopes which varied 
between 5^5 and 100^ of that given by the theory with most 
results close to 705. 



and one -half inch diameter pipes and a frequency range from 
300 to 1200 cycles per second. His results were similar to 
those reported by Rudnick. 

0. B. Wilson, Jr. and David Dies conducted a third and 
similar investigation but used an intermediate size pipe 

q 

^Werth, op.cit. 




*rK 



G, V/erth^ conducted a similar investigation using one 
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(I). 13/16" I.D. ) for the purpose of determining whether the 
tube size influenced the observed rate of attenuation. 
Pressure amplitudes of the order of 0.1 atmosp. eres were 
used with frequencies ranging from I4.O to 180 cycles per 
second. rlo simple dependence of attenuation on tube radius 
was found. £ 2 J 
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5. Description of the equipment 

In order to produce the high intensity sound waves that 
were necessary for the rapid formation of the shock 'waves 
studied, it wrs necessary to cut a stream of compressed air 
with a rotor blade. A compressor capable of delivering up 
to 1^.00 cubic feet of air through a range of pressures up to 
ten pound! per square inch was used to feed air to the 
chopper (Pig 3). The compressor was originally designed for 
and used in testing pressurized cabins of aircraft of the 
U. S. Navy for leaks. In order to meet the fire regulations, 
the gasoline motor was removed and a 15 horsepower electric 
motor was mounted on the frame. This was coupled to the 
shaft of the compressor by a triple v-belt drive using a 
two to ore ratio in order to approximate the rated speed 
for which the compressor was designed. 

The chopper consisted of a plate of ^ inch aluminum 
with four holes cut into it. See Pig. These holes were 
closed and opened by means of a rotor, mounted flush with 
the pl-^te and driven by a 0.5 horsepower motor. When the 
rotating blades of the chopper were in such a position that 
the ports were open, the air was permitted to pass through 
freely producing a condensation of the air in the chamber. 

At a quarter turn later, the holes were closed. The momen- 
tum of the air particles carried them into the tube thus 
producing a rarefaction. It was this series of condensa- 
tions and rarefactions that produced high intensity sound 
waves. 
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Schematic Diagram of the Equipment 




Figure 4 Schematic Diagram of the Rotary Chopper 
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vl / 1 




FIG 6 Siren 
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FIG 7 Air Compressor and Siren before Sound Proofing 
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The w?ves then passed into an adjustable chamber con- 
sisting of a fixed pipe five and one-half inches inside dia- 
meter over which slid another pipe six inches inside dia- 
meter. See Pig 6. An airtight fit was maintained by means 
of two rubber O-rings mounted on the outside of the smaller 
pipe. The proper adjustment of the outer pipe permitted it 
to be placed at a distance eoual to an integral number of 
half wave lengths from the chopper. This assured a maximum 
pressure swing at the entrance to the shock tube. It was in 
this tube that the sound wave was to grow into a shock wave. 

It was originally planned to use a barium titrnate 
transducer mounted on the end of a long probe tube inserted 
into the shock tube from the end, to measure the pressure 
variation of the shock waves. The construction and cali- 
bration of these microphones proved to be a very difficult 
job. After the microphones were constructed, they proved 
to be too insensitive to produce a good trace on the oscillo- 
scope. In order to overcome the weakness of the signal 
produced by the crystal, a system employing a doubly shielded 
conductor and a cathode follower preamplifier was tried. 
Although this increased the response, it was found that the 
system was microphonic and, it was felt, would yield highly 
distorted results. 

The system finally decided upon and used to measure the 
pressure variations was to drill and tap holes in the side of 
the shock tubes at one foot intervals starting four feet from 
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the entrance to the growth tube. An Altec 21-Br-200-3V 
microphone w?s fitted with a l/ 8 " x l>r" adapting probe 
tube and inserted into the holes in succession. The pressure 
was read on a Hewlett-Fnckard I 4 .OO volt-meter. In parallel 
with the volt-meter was an oscilloscope on vhich could be 
viewed the shape of the wave and a counter in order to ac- 
curate ly determine the frequency of the shock wave s. By 
reading the counter before data was taken, it was possible 
to eliminate any error that might have been introduced by 
a variation in the frequencies and a concurrent variation 
in intensity. 

The Hewlett-Packard voltmeter is calibrated in deci- 
bels with respect to one milliwatt at 600 ohms. In order 
to convert this reading to decibels with respect to one 
volt, a correction of 2.2 decibeS had to be subtracted 
from the meter reading (see Appendix A). The meter responds 
to the average rectified voltage and is calibrated in terms 
of the r.m.s. voltage of a sine wave. It has been deter- 
mined that when a saw- tooth is fed into the meter, a cor- 
rection of 11 db. must be added to give the peak-to-peak 
voltage change (See Appendix B). 

The Altec 21-BR-200-3V microphone was calibrated 
against a known standard which in our c^se was an Altec 
21-BR-150 microphone which had been calibrated by the 
manufacturer. The original characteristics were very poor 
because of the large air cushion between the probe and the 
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diaphram of the microphone. When a brass plug mes added 
to the inside of the probe base, the air space was greatly- 
reduced and the characteristics flattened out considerably, 
particularly in the range of frequencies at which the inves- 
tigation was carried out. The calibration curve is shown in 
Figure 9.^ 

The shock tubes were a series of twenty foot seamless 
steel tubes with inside diameters of 0.75 inch, 1-^", l-3/V'» 
2'f" and 2-3/V'. These were each fitted with a plate which 
was bolted to the plate welded on the large growth tube of 
the adjustable chamber. On the far end of the tubes was 
placed an absorption section consisting of a straight length 
of pipe followed by two 90° bends. These sections were 
filled with tapered glass wool in such a manner as to change 
the cross-section of the pipe gradually. In this way, the 
shock wave was rapidly absorbed and reflections were kept to 
a minimum. 
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\ 



Calibration curve of 21-BR-l.^o 
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6. Conduct of the investigation 



This investigation was conducted on two tubes, i.e.; 
the 3/4" I.D. and the 1-3/4" I.D. pipes only because a frac- 
turing of the diaphram of the 21-BR200-3V microphone caused 
suspension of the tests. 

The speed of the chopper was controlled with a variac. 
The microphone was placed in each of the holes in the test 
pipes with the remaining holes plugged with a standard screw. 
Care was taken to prevent the screw from entering the tube 
and thereby distorting the wave form. iReadings were taken 
on the volt-meter only after checking the electronic counter 
to make sure the frequency was constant. 

7. Results of the investigation 
Prom Rudnick's equation [lj 



it can be expected that the attenuation per foot of pipe 
should be frequency dependent. This is borne ^by the graphs 
Pig. 12 and 13. If the equation is rearranged slightly we 
have the expression: 



and a graph of 1/^- vs x/^ should be a linear function with 



14 and lf>, however, the slopes vary with frequency. The 
numbers in parentheses are the value of the ratio of the ob- 
served slope to the theoretical slope. This value we have 
called R. 




± _ j_ 

<f (f. -v 




V j 

a constant slope of The linearity is verified in Pig 
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This indicated that the equation above does not ade- 
quately describe the attenuation. If the refinement used by 



Wilson and Bies 



I . X = 

i & 



Q, * 



(f o 






is used and values for and a£ are substituted, we get 

4-4 s till _ no\/f\>‘ ( x -*.) 4 - 

J zr \-\-j ‘ d ~ K 6^ 

substituting c= /\f 



I 

? 



<r 0 ‘ zr v J i/7 l “a / ^ 

and dividing by (x-x 0 )//\ 



(*-£•) 



r-*- / 
^ r 



no 



i/F 



d 



J 

£a'<- 



_ / 



but ? ~ £ is the observed slope of the graphs on Pig llj. 

> r+> 

and 1 5. If both sides of the equation be divided by ~X~f 



we have on the left side of the equation the expression 
that we have called R. 



R = 1 + 






rv x r 

* fH 



no vA 
*< "oTj? 

Again rearranging _ ^ _*/ 

-- /*,, ~~d~ <■-' 

and plotting a graph of (R - 1) /F vs l/f^ for a given gas 

/7<J 

and pipe radius a straight line with a slope of — — 

.06 a 

should result. The experimental data does not give constant 
values hence one must conclude that the equation used by 
’Wilson and Bies does not adequately describe the situation 
either . 

From Fig 14 and 1 5 we see a definite variation with 
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frequency and comparing, the slopes of similar frequencies 
in the two pipes, it can be seen that the waves attenuate 
faster in the smaller pipe. 

When the ratio of observed to theoretical slopes are 
plotted against frequencies, (Fig 16) it is noted that as 
the frequencies increase, the ratio of observed to theoreti- 
cal slopes decreases to a certain value of the frequency 
after which it rises. This would indicate that some pheno- 
menon has been altered radically. The parallelism of the 
curves after the "critical" frequency is of particular 
interest. 

Further investigation should be conducted to determine 
the shape of this curve more accurately and to determine 
whether this relation holds for other sizes of pipes. 
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APPENDIX A 



SAMPLE CALCULATIONS 

1. BR21-200-3V microphone sensitivity calibrated against 

factory calibrated BR21-150 microphone. 

(a) At 300 cps, BR-21-200 calibrated - 31.5 db below 
BR21-150. 

(b) At 300 cps, BR21-150 was factory calibrated at 
-54 db reference to 1 volt (ms) per dyne per 
cm 2 (rms ) . 

(c) Therefore, the BR21-200 is -5k/ ( -31.5 ) = -85.5 db 
below. 1 volt (rms) per dyne per cm 2- (rms) 

(d) As value in (c) above is rms, -lldb additional 
should be included to get peak to peak values of 
a saw tooth wave,i.e. - 8 5.5 / (- 11 ) = - 95 . 5 db. 
(See appendix b) 

(e) 1 dyne per cm £ = / 7 I; db reference to 0.0002 dynes 
per cm* (= 0 db) 

00 -74 “ 96.5- - 170.5 db reference to 1 volt for 
a Odb Sound Pressure Level (SPL) peak to peak. 
This is the reference number for this frequency. 

(f) As Hewlett-Packard voltmeter used reads 2. 2db 
reference to 1 volt, / 2. 2 db must be subtracted 
from each reading. However, subtracting 2.2db 
from reference 170 . 5 db allows use of readings 
directly. Therefore, corrected reference is 
-170.5 / 2.2 = -168. 3db m on the lewlett-Packard 

Voltme ter, 

* . J- 
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(Note) db r » db's relative to voltage that gives 1 milli- 
watt with 600 ohms resistance. 

2. At station 1 ( 8 1 from siren) for 300 riding cn volt- 
meter » / 9.2 dbj^. 

Therefore, /9.2db - (-168. 3db) = /l77 # 5db m SPL (peak to 
peak,) 

P - P 

3. Using formula, db = 20 log .0002' = 177.5 db . Solving, 

P 2 _p i = 7.41 x 10 8 

~7oM2 

p 2 ~^1 ~ (7»4l x 10 8 ) (2 x 10 ^) s ii (..82 x 10 ^- dynes/cm^ 

1 atmosphere = 1.013 x 10 ^ dynes/cm^ 

6°t> ( p p -P-i ) PK to Pk = 14.82x104 = , 1 ) 4.75 atmospheres 

4. Calculation for 4r , where = ?z. _p i 

* ^1 



also 1 = P]_ 

3 



P, -Pi 



P- -Pi " \ 



‘*‘?r L 



Assuming P Q = 1 atmosphere 

- 2 “ P 1 = PP to pi: pressure of .1475 atmosphere calcu- 
lated by 3 above. 

For station yl , 300 cps, 1" 0D pipe 

} ■ TTT T7F ’ I = 6,?G "* 5 = 6,28 

Note: v of sound at 20°C = 1086' /Sec 
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APPENDIX B 



COMPUTATION FOR CONVERTING RI-'S VALUE OF SAW-TOOTH WAVE INTO 

PEAK TO PEAK VALUE 



1. 1 dyne/cra2 -85.5 db below 1 Volt (for microphone used) 

65.5 = 4.275 

20 



antiloc of 4.275 = 1.88xlo4 
-85.5 db^ l = .53 x io'b 

1.88'(10^)= 53 x- volts 

2. 53/* volts ^ 1 dyne/cm 2 - 

Saw tooth wrve formula, p = ppp 




volts 



(Sin wt-1 sin 2wt/l_ sin3wt 
2 3 




Peak to Peak 
Pressure 



..) 



3. For sine wave 

E » E q sin wt 

avg = 2 E 0 
TT 

•^rms = .° 

\TT 

^rms = 'TT~ ° 1.11 

^avg 2 

-‘-’avg = ^xvns (rerding on voltmeter)^: .92 db 

~T7Tl 

or E aV g 1 db less than E read on Hewlett Packard 



29 



4. Prom 2 above, P avg = 1 P pp 



P PP = h p avg = gges / cn g . E 

53 volts ° 

®avs = E( reading) 

l.lT 

*PP » ^TTl E (reading dynes/cm^ 

P PP = 3.6 E (rending) 
and 3.6^ 11 db 

0 d 0 add 11 db to rms readings to get pe r>k-to-peak 
sound pressure level of saw tooth wave. 
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APPENDIX C 



TABLE SHOWING RATIO OF OBSERVED TO THEORETICAL ATTENUATION 

SLOPES 
1" PIPE 



f 

cps 


* 


Theoretical 

-yt. 


Slope obs. Slope 

for air) (from graph) 


R r Obs 

Theo 


298 


3.65 


.86 


.88 


1.025 


583 


1.86 


.86 


.709 


.827 


730 


1.1*8 


. 86 


.625 


.728 


883 


1.20 


.86 


.61*5 


.751 



ljr" PIPE 



583 


1.86 


.86 


.622 


.722 


710 


1.53 


.86 


.572 


- .665 


900 


1.20 


.86 


.635 


.738 
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